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C^ \ We study fractional configurations in gravity theories and La- 

grange mechanics. The approach is based on Caputo fractional deriva- 
tive which gives zero for actions on constants. We elaborate fractional 
geometric models of physical interactions and we formulate a method 

y\i ' of nonholonomic deformations to other types of fractional derivatives. 

C^ ■ The main result of this paper consists in a proof that for corresponding 

classes of nonholonomic distributions a large class of physical theories 
are modelled as nonholonomic manifolds with constant matrix curva- 
ture. This allows us to encode the fractional dynamics of interactions 
and constraints into the geometry of curve flows and solitonic hierar- 
chies. 
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1 Introduction 

It is well known that both the modified Korteweg - de Vries, mKdV, 
and the sine-Gordon, SG, (solitonic) equations can be encoded as flows of 
the curvatures invariant of plane curves in Euclidean plane geometry [TH El 
El QUI 13 [12] • Such constructions were developed [H [151 Q] f° r curve flows 
in Riemannian manifolds of constant curvature [HI [IS] which gave rise, for 
instance, to a vector generalization of the mKdV equation and encode its 
bi-Hamiltonian structure of multi-component mKdV equations. 

A crucial condition behind above mentioned geometric models was that 
the frame curvature matrix is constant for a series of examples of Riemann 
symmetric spaces, with group symmetries etc, for which certain classes of 
nonlinear/ solitonic equations can be associated. It was considered that the 
approach would have a limited importance, for instance, for study of more 
general classes of curved spaces (pseudo-Riemanian ones, Einstein manifolds 
etc) if such geometries are not with constant curvature matrix coefficients. 
Nevertheless, it was possible to elaborate a similar method of solitonic en- 
coding of geometric data for generalized geometries with non-constant cur- 
vatures using certain ideas and methods from the geometry of nonholonomic 
distributions and modelling Lagrange-Finsler geometries [18] . Such models 
with nontrivial nonlinear connection structure were elaborated on (co) tan- 
gent bundles (see [13] and references therein) and our proposal was to apply 
such methods on (pseudo) Riemannian and Einstein manifolds endowed with 
nonholonomic distributions. 

We developed a geometric formalism when nonholonomic deformations 
of geometric structures determined by a fundamental Lagrange/ Finsler / 
Hamilton generating function (or, for instance, and Einstein metric) induce 
a canonical connection, adapted to a necessary type nonlinear connection 
structure, for which the matrix coefficients of curvature are constant [TSJ 12] ■ 
For such an auxiliary connection, it is possible to define a bi-Hamiltonian 
structure and derive the corresponding solitonic hierarchy. 



As it is known the fractional calculus which deals with derivative and 
integrals of arbitrary orders [221 [231 I2H [2SJ [2JJ] is an emerging field and it 
has many applications in various fields of science and engineering [301 EH 132] 
[331 [231281 Ml M E5]. 

In a series of our works [201 IZH SI E] (we recommend the reader to see 
details, discussions and bibliography, including notation conventions etc), 
we proposed program of geometrization of fractional calculus applications in 
classical and quantum fractional mechanics, field theories and Ricci flow evo- 
lution theories with fractional/noncommutative/ stochastic etc derivatives. 
It should be emphasized here that it is not possible to elaborate an unified 
geometric formalism for all types of existing definitions of fractional deriva- 
tives. Different such operators have their priorities and minuses and result 
in very different geometric (integro-differential) and, for instance, physical 
implications. Nevertheless, there is a class of fractional derivatives resulting 
in zero for actions on constants. This property is crucial for constructing ge- 
ometric models of physical and mechanical theories with fractional calculus 
which (from many formal points of view) are very similar to those on inte- 
ger dimensional spaces but with corresponding nonholonomic distributions 
selecting a prescribed type of fractional dynamics and/or evolution and, in 
general, nonholonomic constraints. Such an analogy allows us to involve in 
our study a number of very powerful geometric methods already used, for in- 
stance, in deformation quantization, geometric mechanics, Lagrange-Finsler 
geometry etc. Via nonholonomic deformations/transforms, the constructions 
based on Caputo fractional calculus seem to have generalizations to include 
other types of fractional derivatives. 

The main goal of this paper is to provide such a geometrization of frac- 
tional calculus formalism [22J[2J[2J[2il|2^[271[2^[2JEa[321E3] extended to 
certain applications/modifications of gravity and Lagrange mechanics theo- 
ries, in a fractional calculus manner, which will allow in the second our 
partner work [6] to encode such geometric constructions into (fractional) bi- 
Hamiltonians and associated solitonic hierarchies. In brief, two of our works 
(this one and [B]) present a fractional version of articles [121 12]. 

The paper is organized as follows: 

In section 2, we outline the geometry of N-adapted fractional spaces. 
In section 3, we provide an introduction into the theory of fractional grav- 
ity, show how the fundamental equations in such theories can be integrated 
in general form and provide a geometrization for fractional Lagrange me- 
chanics. Section 4 is devoted to the geometry of fractional curve flows and 



related models of fractional N-anholonomic Klein spaces. The Appendix 
contains necessary definitions and formulas on Caputo fractional derivatives 
and related nonholonomic fractional differential geometry (with nonlinear 
connections). 

2 Fractional Manifolds with Constant Curva- 
ture Coefficients 

Let us consider a " prime" nonholonomic manifold V is of integer dimen- 

Its fractional extension V is modelled by 

a a a a 

a quadruple (V, N, d, I), where N is a nonholonomic distribution stating a 
nonlinear connection (N-connection) structure (for details, see Appendix [A] 

a 

with explanations for formula (1A.4J) . The fractional differential structure d 
is determined by Caputo fractional derivative (jA.lj) following formulas (1A.2J) 



and (1A.3J) . The non-integer integral structure I is defined by rules of type 



(Q. 



The goal of this section is to prove that it is possible to construct a 
metric compatible linear connection D with constant matrix coefficients of 
curvature, computed with respect to 'N-adapted' frames, for any fractional 

q, a 

metric g (1A.7J) (equivalently, ( 1 A ■ 8 1) ) on a nonholonomic manifold V. 

2.1 N— adapted frame transforms and fractional met- 
rics 

For any respective frame and co-frame (dual) structures, 

a 

a e a > = ( a ei>, a e a i) and "e* 3 = ( a e l , a e a ) on V, we can consider frame 
transforms 

a e a = A a a '(x,y) a e a , and V = A%,(x,y) V. (1) 



1 A nonholonomic manifold is a manifold endowed with a non-integrable (equivalently, 
nonholonomic, or anholonomic) distribution. There are three useful (for our considera- 
tions) examples when 1) V is a (pseudo) Riemannian manifold; 2) V = E(M), or 3) 
V = TM, for a vector, or tangent, bundle on a base manifold M. We also emphasize that 
in this paper we follow the conventions from Refs. [HI [20j [21] when left indices are used 
as labels and right indices may be abstract ones or running certain values. 



A subclass of frame transforms (pQ), for fixed "prime" and "target" frame 
structures, is called N-adapted if such nonholonomic transformations pre- 
serve the splitting defined by a N-connection structure N = {N?}. 

Under (in general, nonholonomic) frame transforms, the metric coefli- 

a a 

cients of any metric structure g on V are re-computed following formulas 
a g a p(x,y) = A a a '(x,y) Af'(x,y) a g a >p>{x,y). 

a a 

For any fixed g and N, there are N-adapted frame transforms when 

| = a gij (x,y) V® a e j + a h ab (x,y) Q e a ® V, 
= a g iff (x,y) V<8> V+ a h a/b ,{x,y) a e a ' ® V, 



where a e a and a e a are elongated following formulas (1A.6J) . respectively by 



"N^ and 



a N a f =A a a (x,y)A^,(x,y) a N a 3 (x, y), (2) 

or, inversely, 

a N a 3 = A al a (x,y)A^y) a N a ' f (x,y) 

with prescribed a N a -,. 

We preserve the N-connection splitting for any frame transform of type 
(P) when 

Oi _ A i A3 o. 01 A a A b a h 

yi'j' — -n- i'Sl j, yij } ll'a'b' — -n- a ;/1 b' ' l abi 

for A i % constrained to get holonomic a e l = A i % a e\ i.e. [ a e % , a e J ] = and 
a e a = dy a + a N a -idx^ , for certain x l = x l (x\ y a ) and y a = y a (x\ y a ), with 
a N a -, computed following formulas (T5]). Such conditions can be satisfied by 
prescribing from the very beginning a nonholonomic distribution of necessary 
type. The constructions can be equivalently inverted, when a g a i3 and a Nf 
are computed from a g a 'i3> and a N°t , if both the metric and N-connection 

a 

splitting structures are fixed on V. 

2.2 d— connections with constant curvature coefficients 

From the class of metric compatible fractional d-connectionsj uniquely 

ex. 

defined by a fractional metric structure g, we chose such a N-connection 



2 see definitions and main formulas in Appendix lAl 



splitting with nontrivial coefficients a N^(x, y) when with respect to a N- 
adapted frame the canonical d-connection (IA.140 has constant coefficients. 

a a 

Any fractional metric g on V defines a set of metric compatible fractional 
d-connections of type 

a f 7 ' - ( a p' -n a f a ' -- a T a ' - rrm<it a r i ' — n *r°' — r\\ 

1 a' /3' ~ I U 3'k' — u ) u Vk' — O^b'k' ~ COTlbl, ^>j, c , — U, O b / c / — U I 

(3) 
with respect to correspondingly constructed N-adapted frames (1A.5I) and 
(1A.6J) . when a N = { a Nf, (x,y)} is a nontrivial solution of the system of 
equations 

2 il$ v = %{ <*N$) - a h a ' c ' «h d , v h a N d k , (4) 

for any nondegenerate constant-coefficients symmetric matrix a hd'b' an d its 

inverse a h a c . The Caputo fractional derivative d b , allows us to dub such 
constructions for fractional spaces [ISJ |2| . 

For both integer and non-integer dimensions, the coefficients "r 7 ^,™ of 

the corresponding to g Levi-Civita connection "V are not constant with 
respect to N-adapted frames^ 

The curvature tensor of fractional d-connection q r 7 Q//3 , (j3j) defined by a 

metric g has constant coefficients with respect to N-adapted frames a e a / = 
[ a e if , a e a >] and a e a ' = [ V, a e a '|^with a Nf, subjected to conditions (0}. 
Introducing constant coefficients oT 7 ^, © into formulas (1A.10D . we get 

art a' /ap*' n ct pa' arc' a j a' arc' a r a! 

JX p, 7 , S , — { o-rt h ij, k i — U, o-ft b ,ji k , — q±j V ji qIj c i k , - iv b i k i qL, c ,j, 

= COnst, P h'j'a' = 0; qP b'j'a 1 = ®i $ j'b'c' = 0) O 1 ^ b'd'c' = 0)- 

Such proprieties hold for a class of prescribed nonholonomic distributions. 
Of course, in general, with respect to local coordinate (or other N-adapted) 
frames, the curvature d-tensor a H a a s does not have constant coefficients. 
Nevertheless, we can always chose certain nonholonomic N-connection config- 
urations when the corresponding adapting will generate d-connections with 
constant matrix curvature coefficientso 



3 In explicit form, such coefficients are computed following formulas (|A.f5|) . 
Using deformation relation (|A.15j) . we can compute the corresponding Ricci tensor 
"R°0 s for the Levi-Civita connection "V, which is a general one with 'non-constant' 



coefficients with respect to any local frames. 



Using formulas ( lA.lip and (IA.12J) . we compute 

a W _j_ a rr a l fi' art o.J/j' an i aua'b' an 

**■ — § K-a'P — 9 % + ft 'Vfc' 

= ^+ ^ = const. 

We conclude that a fractional d-connection qT 1 ^^ (J3J) has constant scalar 
curvature. 



3 Fractional Einstein and/or Lagrange Spaces 

The goal of this section is to show how the integer dimensional Einstein 
gravity and Lagrange mechanics can be encoded and/or generalized in terms 
of fundamental geometric objects on nonholonomic fractional manifolds. We 
review and develop our constructions from J20J, [211 HI E] • 

3.1 Fractional gravity 

An unified approach to Einstein-Lagrange/Finsler gravity for arbitrary 
integer and non-integer dimensions is possible for the fractional canonical d- 
connection a D. The fractional gravitational field equations are formulated 
for the Einstein d-tensor (1A.13[) . following the same principle of constructing 
the matter source "Y^ as in general relativity but for fractional metrics 
and d-connections, 

a E ps = a Y^5. (5) 

Such a system of integro-differential equations for generalized connections 
can be restricted to fractional nonholonomic configurations for a V if we 
impose the additional constraints 

a L c a] = a e a ( a N<) } a d} b = } «n» ii = 0. (6) 

There are not theoretical or experimental evidences that for fractional 
dimensions we must impose conditions of type (Q but they have certain 
physical motivation if we develop models which in integer limits result in the 
general relativity theory. 



3.1.1 Separation of equations for fractional and integer dimen- 
sions 



We studied in detail [21] what type of conditions must satisfy the coeffi- 
cients of a metric (1A.8J) for generating exact solutions of the fractional Ein- 
stein equations (jSJ). For simplicity, we can use a "prime" dimension splitting 
of type 2 + 2 when coordinated are labelled in the form u 13 = (x\y 3 = v, y 4 ), 
for i,j, ... = 1, 2. and the metric ansatz has one Killing symmetry when the 
coefficients do not depend explicitly on variable y 4 . 

The solutions of equations can be constructed for a general source of typqj 



t T a l3 = diag[ a T y ; a T 1= a T s 



"T 2 (x k ,v) 



*r ? 



aj 



*r 4 (x h )} 



13 — """*yi '71 

For such sources and ansatz with Killing symmetries for metrics, the Einstein 
equations (J5J) transform into a system of partial differential equations with 
separation of equations which can be integrated in general form, 

x\ a g? 9l 92 



a ol 



Ri 



( P 2 

Il9 



2 a gi a g 2 



"92 



"92 



+ V 



a d3 
K 3 



a p4 



a K) 



2 a h 3 a h 4 l 
a h* 3 a hl, 



a g[ % 

2 a 92 

\ a hT 



2 a 9i 
2 a 9i 



Ey 



4, 



2 a h d 



2 «h« 



J-2, 



(7) 



(8) 



"R 



.u 



2 a h d 



x h\ 



a h* A 



2 a h A 



2 a h 3 



(9) 



a K 



lX id_ x i a h 3 dj^ 



lX kdj.k a h\ 



*h* 



hi 



a R. 



U: 



a h A 



n l 



+ 



2 a h 3 K \ a h 3 
where the partial derivatives are 



"*4 ai* " ai* 

'h 9 4 



(10) 



!x 2 d x 2 a a, 



d v a = 1V d v a a, 



a a* = d x a = ix i d x i a a, a a' = d 2 a 

°such parametrizations of energy— momentum tensors are quite general ones for various 
types of matter sources 



8 



being used the left Caputo fractional derivatives (1A.2J) . 

Configurations with fractional Levi-Civita connection Q V, of type 
can be extracted by imposing additional constraints 



w* = a eM\ a h±l a e k a w i= a e, a w k , 

< = 0, a d l a n k = d. k a n l . (11) 



c„* 



We can construct 'non-Killing' solutions depending on all coordinates 
when 

a g = a g t (x k ) a dx i ® a dx i + a u 2 (x j ,v,y 4 ) a h a (x k ,v) a e a ® a e a , 
V = a dy 3 + a w t (x\v) a dx\ a e 4 = a dy 4 + a n t (x k , v) a dx\ (12) 

for any a w for which 

a e k a u = d k a u+ a w k a u* + a n k d y4 a u = 0. 

3.1.2 Solutions with a h* 3 4 ^ and a T 2 , 4 ^ 

For simplicity, in this paper we show how to construct exact solution 
with metrics of type f fl2|) when a h 34 ^ (in Ref. [21], there are analyzed 
all possibilities for coefficientqj) We consider the ansatz 

a g = e aip{xk) a dx i ® a dx i + h 3 {x k , v) Q e 3 ® °e 3 + h 4 {x k , v) a e 4 ® a e 4 , 
Q e 3 = a dv+ a Wi{x k ,v) a dx\ a e 4 = a dy 4 + a n i {x k ,v) a dx i (13) 

when 

<V+ >" = 2 Q T 4 (x fc ), (14) 

a h\ = 2 a h 3 a h 4 a T 2 {x i ,v)/ a <f)*, (15) 



a /3 a Wi + a ai = 0, (16) 

X* + a 7 X = 0, (17) 

where > = In I , *M =1, a 7 = (in \ a h 4 \ 3/2 /\ a h z \)* , (18) 
a ai = a h{d k a 0, a (3 = a h\ Q 0* . 



i; 



by nonholonomic transforms, various classes of solutions can be transformed from one 



to another 



For a h\ ^ 0; a T 2 ^ 0, we have a (p* ^ 0. The exponent e a ^ xk) is the frac- 
tional analog of the "integer" exponential functions and called the Mittag- 

a 

Leffler function E a [(x— 1 x) a }. For a ip(x) = E a [(x— 1 x) a ], we have d_iE a = E a , 
see (for instance) [17] . 

Choosing any nonconstant a <p = a <p(x l ,v) as a generating function, we 
can construct exact solutions of ffT4l) - ffT7|) . We have to solve respectively the 
two dimensional fractional Laplace equation, for a gi = a g 2 = e a ^ x >. 
Then we integrate on v, in order to determine "h^, a h 4 and a rii, and solving 
algebraic equations, for a Wi. We obtain (computing consequently for a chosen 
a <f)(x k ,v)) 

a 9i = a 92 = e "«**>, «% = ± i^^i, (19) 



"To 



X W; = -d t a 6l a J>\ 



a n t = a ink (x l )+ > fc (0 lV Iv[ a h 3 /(y/\^\)% 

a 

where Qh^(x k ), "n k (x l ) and 2 n k {x l ) are integration functions, and lV I v 
is the fractional integral on variables v. 

We have to constrain the coefficients (fT9l) to satisfy the conditions ( ITT]) 
in order to construct exact solutions for the Levi-Civita connection a V. 
To select such classes of solutions, we can fix a nonholonomic distribution 
when 2 n k (% l ) — and 1n k (x l ) are any functions satisfying the conditions 

a a 

d_ t fnk (xj) = d k iiii (x J- ) . The constraints on a (fi(x k ,v) are related to the 

a 

N-connection coefficients a Wi = —d t a <p/ a <p* following relations 



( a w t [ >])* + « Wi [ >] ( a h 4 [ >])* + d t a K[ >] = 0, 

d i a w k [ a (f ) ] = d k a Wi [ a <j>l (20) 

where, for instance, we denoted by a hi[ a 4>] the functional dependence on 
a <p. Such conditions are always satisfied for a cf) = a <j)(v) or if a <j) = const 

when a Wi(x k , v) can be any functions as follows from (TT6]) with zero a f3 and 
a a i) see ([T~ 
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3.2 A geometric model for fractional Lagrange spaces 

Any solution a g = { a g a 'p'(u a )} of fractional Einstein equations (j5]) can 
be parametrized in a form derived in previous section. Using frame trans- 
forms of type a e a = e a ' a a e a <, with a g a/3 = e a ' a e^ a g a '/3>, for any a g a(3 (|A.8p 



we relate the class of such solutions, for instance, to the family of metrics of 
type (fT31) . Such solutions can be related also to analogous fractional models, 
via corresponding e a a to a Sasaki type metric, see below ( ]2T]) . 

a 

A fractional Lagrange space is defined on a fractional tangent bundle T_M 

a a. a 

of fractional dimension a 6 (0,1) by a couple LF_ = (M,L), for a regular 

a a 

real function L : T_M — > M, when the fractional Hessian is 



I / a a a a \ a 

l 9ij = ^ ( 3^ + 3$ J L ^ 0. 



Any Z^_ can be associated to a prime "integer" Lagrange space L n . 

Let us consider values y k (r) = dx k (r)/dT, for x(r) parametrizing smooth 
curves on a manifold M with r e [0, 1]. The fractional analogs of such con- 
figurations are determined by changing d/dr into the fractional Caputo 

q a a a 

derivative d T = lT <9 r when a y k (r) = d T x k (r). Any L defines the fundamen- 
tal geometric objects of Lagrange spaces when 

a a a a a 

1. the fractional Euler-Lagrange equations d T ( ^yid^L) — lX i d_ t L = are 
equivalent to the fractional "nonlinear geodesic" (equivalently, semi- 



2 



spray) equations I d T x + 2G (x, a y) = 0, where 



a i a 

G k = - L g k > 



a a \ a a 

V \yiQ-j I ix l <liL 1 - 1 x i <liL 



defines the canonical N-connection "iV? = iy jdjG (x, a y); 
2. the canonical (Sasaki type) metric structure is 

lI = a L g k j(x, y) a e k ® V + a L g cb (x, y) £e c ® £e fe , (21) 

where the frame structure 2 e ^ = ( 2 e «) e a) is linear on 2-/V"; 

11 



3. the canonical metrical d-connection 

«D = (h «D,v «D) = { «T\p = ( a V jk , «C) C )} 

is a metric compatible, "D I ^g ) = 0, for 

ayii ayii Q ~7 T i k _i_ /^i a ^ c 

c L j — c L j-y L e — n jk e I" °jc L e 5 

with D jk = L\ k , % = CI in «T\ = «TX tef = L\ k e k + Cfc £e c Q 
and generalized Christoffel indices 

«n _ o Jvr I a a , a a a a \ 

n jk — 9 L ^ \ L k L ^ r L ■? L ^ fcr " L r Lyjk) ) 

<^bc - 2 l9 { e c L 5f M + e c L # cd - e d L g bc ) . 

4 Nonholonomic Fractional Curve Flows 

We formulate a model of geometry of curve flows adapted to a N-connection 

a 

structure on a fractional manifold V. 

4.1 Non— stretching and N— adapted fractional curve flows 

A canonical d-connection operator a D ( IA.14J) acts in the form 

a D x a e a = (XJ Q rj) a e 7 and a B Y a e a = (YJ a TJ) a e 7 , (22) 

where " |"' denotes the interior product. The value a Dx = X' 3 "D^ is a 
covariant derivation operator along curve j(t, 1). It is convenient to fix the 
N-adapted frame to be parallel to curve 7(1) adapted in the form 

V = hX, for i = 1, and V, where h a g(hX, V) =0, (23) 

a e n+1 = vX, for a = n + 1, and a e s , where v a g(vX, a e s ) =0, 

for i — 2, 3, ...n and a = n + 2, n + 3, ..., n + m. The covariant derivative of 
each "normal" d-vectors a e a results into d-vectors adapted to 7(7", 1), 

Q D X V = -p 7 (u) X and a D hX .hX = p\u) a ej, (24) 

a D x Q e s = -p a (w) X and a B v ^vX = p s (w) Q e a , 



7 for integer dimensions, we contract "horizontal" and "vertical" indices following the 
rule: i = 1 is a = n + 1; i = 2 is a = n + 2; ... i = n is a = n + n" 

12 



which holds for certain classes of functions p % {u) and p a (u). The formulas 
( T2"2"j) and (1241) are distinguished into h- and v-components for X =hX. + 
vX and Q f>= {KD,vD) for a D = { a f 7 Q/3 },/iD = { a L) kl a L a bk ] and 

vB = { «dj e , a cs c }. 

a 

A non-stretching curve 7(7", 1) on V, where r is a real parameter and 1 

a 

is the arclength of the curve on V, is defined with such evolution d-vector 
Y = 7 r and tangent d-vector X = 71 that a g(X, X) =1. Such a curve 

a a 

7(r, 1) swept out a two-dimensional surface in T/nV c TV. Along 7(1), 
we can move differential forms in a parallel N-adapted form. For instance, 
°T x =f XJ a T a P. Such fractional spaces can be characterized algebraically 
if we perform a frame transform preserving the decomposition (1A.4J) to an 
orthonormalized basis a e a i, when 



a e„ ->■ e ° (u) 



(25) 



is an orthonormal d-basis. In this case, the coefficients of the d-metric (1A.8J) 
transform into the (pseudo) Euclidean one a g a /^/ = r) a >p' by encoding the 
fractional configuration into the structure of local d-bases. We obtain two 
skew matrices 



T& = hX\ °T 



a-pt 3 



e h*P J 



and 



aj-\a b 
1 uX 



aria/) 



wxj a r 



2 e 



^xP 



^ x = g(hX,e l ) = [l } 1 _^ and e^ = h(vX,e a ) = il,0 1 _ : ^0) 1 

n— 1 m— 1 



- 1 /iXi' 



p> 
-Pv 0\h 



and "f\ Xc / 



p v 

—Pa' 0[v] 



with 0[h] and 0^] being respectively (n 
matricesjj 



1) x (n — 1) and (ra — 1) x (ra — 1) 



8 The above presented row-matrices and skew-matrices show that locally an N- 

anholonomic fractiona manifold V are related to prime spaces of integer-dimension n + m. 
With respect to distinguished orthonormalized frames the constrained fractional dynamics 
is characterized algebraically by couples of unit vectors in R™ and M. m preserved respec- 
tively by the SO(n — 1) and SO(m — 1) rotation subgroups of the local N-adapted frame 



structure group SO(n) SO(m). The connection matrices 



A hX i> 



and 



uX< 



be- 



long to the orthogonal complements of the corresponding Lie subalgebras and algebras, 
so(n — 1) C so(n) and so(m — 1) C so(m). 
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The torsion and curvature tensors (lA.lOh can be written in orthonor- 
malized component form with respect to ( 12 3 p mapped into a distinguished 
orthonormalized dual frame ( 125]) . 

a fa' ^ ag x ^ _ afy^ a^ + a^ af^a> _ a^ a f^of ( 2 6) 

and Q 7^'(X,Y) = a D Y a f x / - Q D X Q f Y / 

i Q-p 7 a-p a' Q-p 7 «p a' (07^ 

"+" - 1 Y/3' - 1 X 7 ' ~~ x X/9' - 1 Y 7 ' • l z ' J 

The values Q e Y ' = Q g(Y, a e a '), a f Y / = YJ a f ^ = a g( Q e a ', a D Y a e^) 
define respectively the N-adapted orthonormalized frame row-matrix and 
the canonical d-connection skew-matrix in the flow directs, and 



'■■ 



a 1l$ (X, Y) = a g( a e a , [ ft D x , a D Y ] a e pi ) is the curvature matrix. 

4.2 N— anholonomic fractional manifolds with constant 
matrix curvature 

The geometry of integer dimensional Einstein and Lagrange-Finsler spa- 
ces can be encoded into bi-Hamilton structures and associated solitonic hi- 
erarichies [HI [2]. The goal of this section is to show that there is a geo- 
metric background for extending the constructions to the case of fractional 
spaces. We shall elaborate the concept of fractional N-anholonomic Klein 
space which in [6] will be applied for constructing fractional solitonic hierar- 
chies. 

4.2.1 Symmetric fractional nonholonomic manifolds 

For trivial N-connection curvature and torsion but constant matrix curva- 
ture on spaces of integer dimension, we get a holonomic Riemannian manifold 
and the equations (1261) and (I2~T|) directly encode a bi-Hamiltonian structure 
[T5l [TJ. A well known class of Riemannian manifolds for which the frame 
curvature matrix constant consists of the symmetric spaces M = G/H for 
compact semisimple Lie groups G D H. A complete classification and sum- 
mary of main results on such integer dimension spaces are given in Ref . [9] . 
Using the Caputo partial derivative, such a classification can be provided for 
fractional spaces of constant matrix curvature. The derived algebraic classifi- 
cation is that for a used "prime" integer dimension space but the differential 
and integral calculus are those for fractional nonholonomic distributions. 
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Algebraically, our aim is to solder in a canonic way the horizontal and 
vertical symmetric Riemannian spaces of dimension n and m with a (total) 
symmetric Riemannian space V of dimension n + m, when V = G/SO(n + m) 
with the isotropy group H = SO (n + m) D 0{n + m) and G = SO{n + m+\). 
The Caputo fractional derivative is encoded for constructing tangent spaces 
and related geometric objects. For the just mentioned horizontal, vertical 
and total symmetric Riemannian spaces one exists natural settings to Klein 
geometry. A prime fractional metric tensor h a g = { a g%j} on h "V is defined 
by the Cartan-Killing inner product < •, • >h on T x hG ~ hg restricted to the 
Lie algebra quotient spaces hp =hgjh\), with T x hH ~ hi), where hg =hi)®hp 
is stated such that there is an involutive automorphism of hG under hH is 
fixed, i.e. [/i(),/ip] C hp and [/ip,/ip] C hi). In a similar form, we can define 
the group spaces and related inner products and Lie algebras, 

for vg = {h ab }, < •, • >„, T y vG ~ vg, vp =vg/vi), with 
T_ y vH ~ vt),vg =t>f) ©fp,where [ff),fp] Q vp, [vp,vp] C vty, 

for g = {g a p}, < -, • > g , Z {x , y) G ~ g, p = g/i), with 
T {Xty) H ~ fj,fl = f)©p,where [\), p] C p, [p, p] C f). 

Similar formulas in [191 E] are for usual partial derivatives with not underlined 
symbols T. We parametrize the metric structure with constant coefficients 
on V = G/SO(n + m) and fractional differentials in the form 

a g= a g^(dvX) a ®(du^-) a , 



where the coefficients of a fractional metric of type (1A.7J) are parametrized 
in the form 



a ° 

9af3 



% + Q iVf Q A^ a h ab a N* a h ae 



a N[ a h be a h 



ab 



with trivial, constant, N-connection coefficients computed a A e = a h eb a gj b 

for any given sets h eh and a (jj b , i.e. from the inverse metrics coefficients 
defined respectively on hG = SO(n + 1). We can also define an equivalent 
d-metric structure of type (1A.8I) 

°g = a ° gij a e <g) « e j + a A a6 a e a ® a e b , 
e l = (dx l ) a , a e a = (dy a ) a + "N^dx 1 ) . 
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We note that for integer dimensions such trivial parametrizations de- 
fine algebraic classifications of symmetric Riemannian spaces of dimension 
n + m with constant matrix curvature admitting splitting (by certain al- 
gebraic constraints) into symmetric Riemannian subspaces of dimension n 
and m, also both with constant matrix curvature and introducing the con- 
cept of N-anholonomic Riemannian space of type a V = [hG = SO(n + 1), 
vG = SO{m + 1), a Nf\. Such spaces of constant distinguished curvature 
are constructed as trivially N-anholonomic group spaces which possess a Lie 
d-algebra symmetry so^(n + m) = so(n) © so(m). 

A fractional generalization of constructions is to consider nonholonomic 
distributions on V = G/SO(n+m) defined locally by arbitrary N-connection 
coefficients a N°-(x,y) with nonvanishing a W^a and a fl^ but with constant 
d-metric coefficients when the fractional metric Q g is of type a g a >p> = 
[ a 9i'j'i a h a 'b'] (IA.8J) with constant coefficients a gvji = ogi>j> = a g%j and 



a h a 'b' = o^a'b' = a h a b (in this section induced by the corresponding Lie 
d-algebra structure S0ft(n + m)). Such spaces transform into N-anholonomic 
fractional manifolds Q V N = [hG = SO(n + 1), vG = SO(m + 1), a N[] 
with nontrivial N-connection curvature and induced d-torsion coefficients of 
the canonical d-connection_|. 

4.2.2 Fractional N— anholonomic Klein spaces 

We can characterize curve flows (both in integer and fractional dimen- 
sions) by two Hamiltonian variables given by the principal normals h v 
and V, respectively, in the horizontal and vertical subspaces, defined by 
the canonical d-connection Q D = (hD,vD), h v = "D/jx/jX = v % Q e^ 
and v v = Q D„x^X = v a a e^, see formulas 023 j) and ( |24p . This normal 
fractional d-vector Q v = ( h z/, v u), with components of type a v a = (z/, 
v a ) = (z/ 1 , u l , z/ n+1 , z/ a ), encoding the Caputo fractional derivative is in the 
tangent direction of curve 7. It can be also considered the principal nor- 
mal d-vector a w = ( h w, v w) with components of type a w a = (ro 8 , 
w a ) = {xu , vo'\ vj n , w a ) in the flow direction, with 



9 see formulas (|A.10j) computed for constant d-metric coefficients and the canonical 
d-connection coefficients in (|A.f4j) 
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representing a fractional Hamiltonian d-covector field. We argue that the 
normal part of the flow d-vector hj_ = Yj_ = h l a e^ + h a a e-& represents 
a fractional Hamiltonian d-vector field and use parallel N-adapted frames 
a e a / = ( a eii, a e a i) when the h-variables z/', w % \ h 1 ' are respectively encoded 
in the top row of the horizontal canonical d-connection matrices °T hXi , J 
and Q r hYi , 3 and in the row matrix (ey) ± 4= e l Y — g\\ e l x , where g\\ = 
a g(hY,h~X.) is the tangential h-part of the fractional flow d-vector o 

A N-connection structure (in particular, a fractional Lagrangian) induces 
a N-anholonomic Klein space stated by two left-invariant hg- and v g-valued 
Maurer-Cartan form on the Lie d-group G = (hG, vG) is identified with the 
zero-curvature canonical d-connection 1-form qT = { qPo,}, 

a -pa' a -pa' 7' a j i' a k' 1 a /~ii' a k' 

G L /3' — G L fi'yV — hG-^ j'k' e "r uG b j'fe' e • 

For n = m, and canonical d-objects (N-connection, d-metric, d-connection, 
...) derived from (IA.8I) . and any N-anholonomic space with constant d- 
curvatures, the Cartan d-connection transform just in the canonical d-connection 
(1A.14J) . Using the Lie d-algebra decompositions g = hgQ)vQ, for the horizon- 
tal splitting: hg = so(n) © hp, when [hp,hp] C so(n) and [so(n),hp] C hp; 
for the vertical splitting vg = so(m) © fp, when [wp,t>p] C so(m) and 
[so(m),up] C vp, the Cartan d-connection determines an N-anholonomic 
Riemannian structure on the nonholonomic bundle a E = [hG = SO(n + 1), 
vG = SO(m + 1), Nf\. It is possible to consider a quotient space with 
distinguished structure group "Vn = G/SO(n)@ SO(m) regarding G as 
a principal (SO(n) © SO (m))-bundle over a E, which is a N-anholonomic 
bundle. In this case, we can always fix a local section of this bundle and 
pull-back qT to give a (hg © v g)-valued 1-form "f in a point u E a E. [^J 



10 It is possible to encode v- variables v a , vj a ' , h a in the top row of the vertical canonical 
d-connection matrices Q I\ Xa , 6 and a F vYa i b and in the row matrix (eYj == e^ — 

/i|| e x where h\\ = a h(vY,vX.) is the tangential v-part of the flow d-vector. In a compact 
form of notations, we shall write v Q and w a where the primed small Greek indices 
a',/3', ... will denote both N-adapted and then orthonormalized components of geometric 
objects (d-vectors, d-covectors, d-tensors, d-groups, d-algebras, d-matrices) admitting 
further decompositions into h- and v-components. 

11 There are involutive automorphisms ha = ±1 and va = ±1, respectively, of hQ 
and vq, defined that so(n) (or so(m)) is eigenspace ha = +1 (or va = +1) and hp 
(or vp) is eigenspace her = — 1 (or va = —1). We construct a N-adapted fractional 
decomposition taking into account the existing eigenspaces, when the symmetric parts 
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It is possible to generate a G( = hG © t>G)-invariant fractional d- 
derivative a D with restriction to the tangent space T "Vn for any N- 
anholonomic curve flow j(t, 1) in "Vn = G/SO(n)(B SO(m) is 

a D x Q e = [ a e, 7l J a T] and a D Y Q e = [ a e, 7r J a T] , 

admitting further h- and v-decompositions. The derivatives a D x and a Dy 
are equivalent to f )22|) and obey the Cartan structure equations fl26|) and fl27|) . 
We consider a N-adapted orthonormalized coframe a e a = ( a e J , a e a ) 
identified with the (hp © t> p)-valued coframe a e in a fixed ortho normal basis 
for p =hp © vp Chg © t>g. For the kernel/ cokernel of Lie algebra multipli- 
cations in the h- and v-subspaces, respectively, [ Q e ftX , -} h and [ a e„ x , •]„„ , 
we can decompose the coframes into parallel and perpendicular parts with 
respect to a e x , a e = (e c = he c + ve c ,e c ± = he c ± + ve c ±), for p( = 
hp © wp)-valued mutually orthogonal d-vectors e^ and e c ±, when there 
are satisfied the conditions [e x ,ec] B = but [e x ,e c x] 7^ 0; such conditions 
can be stated in h- and v-component form, respectively, [/ie x , hec] hs = 0, 
[hex,he c ±] hs ^ and [ve x ,ve c ] VB = 0, [ve x , ve c ±] vg 7^ 0. There are decom- 
positions 

T U V N ~ p =hp © vp = q =hg © vg/so(n) © so(m) 
and p =p c ©p c x = {hpc ® vpc) ® {hp c ± ®vp c ±), 

with p|| C p c and p c ± C p ± , where [p||,pc] = 0, < p c ±,p(7 >= 0, but 
[p||, pcA ¥" (i- e - Pc is the centralizer of e x in p =hp © vp C/ig © v g); in h- 
and v-components, one have hp\\ C hpc and hp c ± C /ipj_, where [/ip||, /ipc] — 
0, < hp c ±,hpc >= 0, but [/ipy, /ipc±] 7^ (i.e. /ipc is the centralizer of 
e/ lX in hp Chg) and fpy C vpc and fpc^ <= ^P_i_, where [fp||,fpc] = 0, 
< vp c ±,vpc >= 0, but [vp||,up c x] 7^ (i.e. vpc is the centralizer of e„ x in 
vp Cvg). 



"f =i ( %T+a ( gf)) , with respective h- and v-splitting Q L=| ( £ B L+/io- ( £ B L)) and 

"C=^( " fl C + ha( " fl C)). This defines a (so(n) 0so(m))-valued d-connection fractional 
1-form. The antisymmetric part e =| ( e T—a( e T)) , with respective h- and v-splitting 
h °e =| ( f iS e-h(T ( g fl e)) and w Q e =|( ^e - /ur( ™ fl e)), defines a (hp © vp)-valued N- 
adapted coframe for the Cartan-Killing inner product < ■, ■ > p on T_ U G ~ hg(Bvg restricted 
to T u Vn ~ p. This inner product, distinguished into h- and v-components, provides a 
d-metric structure of type "g = [ a g, a h] (|A.8|) .where Q g =< h a e&h a e >hp and 
a h =< v a e®v a e > vp on Q V N = G/SO(n)® SO{m). 
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Acting with the canonical d-connection derivative Q Dx of a d-covector 
perpendicular (or parallel) to "ex, we get a new d-vector which is paral- 
lel (or perpendicular) to Q ex, i.e. "Dec G pc± (or a ~De c ± G pc). In 
h- and v-components, such formulas are written a 'Df l ^hec G hp c ± (or 
&hxhe c ± G hpc) and "D^x^e^ G vp c ± (or "D^x^c 1 - G fpc)- All such 
d-algebraic relations can be written in N-anholonomic manifolds and canon- 
ical d-connection settings, for instance, using certain relations of type 

a D x (e a ') c = /;,(eO c i and Q D x (e"V = -v^(e^') c , 

for some antisymmetric d-tensors v a ^ = — v' 3 a . This is a N-adapted (SO(n) © SO{m))- 
parallel frame defining a generalization of the concept of parallel frame on 
N-adapted fractional manifolds whenever pc is larger than py. If we sub- 
stitute a e a = ( a e l , a e a ) into the last formulas and considering h- and 
v-components, we construct 5'0(n)-parallel and SO (m)-parallel frames. 

5 Conclusions 

The use of the fractional calculus techniques in differential geometry is 
still at the beginning of its application. The fractional operators reveal a 
complex structure and possess less properties that the classical ones. From 
these reasons the fractional operators are not easy to be used within the 
differential geometry and its applications. In this manuscript the calculations 
were done within Caputo derivative. It was shown that for corresponding 
classes of nonholonomic distributions a large class of physical theories are 
modelled as nonholonomic manifolds possessing constant matrix curvature. 
As a result we encoded the fractional dynamics of interactions and constraints 
into the geometry of curve flows and solitonic hierarchies. 

A Fractional Caputo N— anholonomic Mani- 
folds 

For spaces of fractional dimension, it is possible to construct models of 
fractional differential geometry similarly to certain corresponding integer di- 
mension geometries if the Caputo fractional derivative is used [201 EI] • Ap- 
plying nonholonomic deformations, the constructions can be generalized for 
another types of fractional derivatives. 
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A.l Caputo fractional derivatives, local (co) bases and 
integration 

The fractional left, respectively, right Caputo derivatives are defined by 
formulas 

i/W := pT^y J(x- aO- " 1 (^7) S /(*W (A.l) 



!■'■ 



^^M-n^y/^-^-'f-^l/M^. 



We can introduce d := (dx 3 ) a odj for the fractional absolute differential, 
where dx 3 = (dx 3 ) a ^J_, if \x % = 0. Such formulas allow us to elaborate 

a 

the concept of fractional tangent bundle TM, for a G (0, 1), associated to a 
manifold M of necessary smooth class and integer dimM = n\ 12 \ 

Let us denote by L z ( \X, 2 x) the set of those Lesbegue measurable func- 

tions / on [ jx, 2 #] when ||/|| 2 = (J \f{x)\ z dx) l ^ z < 00 and C z [ ix, ?,x\ 

be the space of functions which are z times continuously differentiable on 
this interval. For any real-valued function f(x) defined on a closed interval 

a 

[ ix, 2 x], there is a function F(x) = lX I x f(x) defined by the fractional 

a x 

Riemann-Liouville integral lX I x f(x) := pAr f(x — x') a ~ l f(x')dx', when 

a 

f(x) = lX d x F(x), for all x G [ ix, 2 x], satisfies the conditions 

a 



!xd x ( lX I x f(x)j = f(x), a > 0, 
J x ( lX d x F(x)) = F(x) - F( ix), < a < 1. 



We can consider fractional (co) frame bases on TM. For instance, a frac- 
tional frame basis e^ = e 'rJvP)d_a, is connected via a vierlbein transform 



12 For simplicity, we may write both the integer and fractional local coordinates in the 
form u 13 = (xi ,y a ). We underlined the symbol T in order to emphasize that we shall 
associate the approach to a fractional Caputo derivative. 
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e Ju@) with a fractional local coordinate basis 



o 



dp = (&]' = ix*' &j'i$v = iy b ' @v) , (A.2) 

for f = 1, 2, ..., n and 6' = n + 1, n + 2, ..., n + n. The fractional co-bases are 
related via e = ei, (u^)du^ , where 



"dvF = (W') a , (dy a ') Q ) . (A.3) 



The fractional absolute differential d is written in the form 



a a . . i x . 



,j\l-a 



d := (dx 3 ) a d j: where dx 3 = (dx- 



3\a_ 



-' ' r(2-a) 



where we consider \x % = 0. The differentials dx 3 = (dx 3 ) a x are used as local 
coordinate co-bases/-frames or the "integer" calculus. For < a < 1, we 

a 

have dx = (dx) l ~ a (dx) a . The "fractional" symbol (dx 3 ) a is related to dx 3 and 
can be used instead of " integer" dx 1 for elaborating a co-vector/differential 
form calculus. Following the above system of notation, the exterior fractional 
differential is 

n 

d = J2 r ( 2 - a)(zT _1 dx j dj. 

3=1 

The fractional integration for differential forms on interval L — [ ix, 2^] 
is performed following formula 



L I[x] lX d x f(x) = /( 2 x) - /( ix), 



when the fractional differential of a function f(x) is lX d x f(x) = [...], when 

2X 

{dx) ~ a y^\(dxT JwOOl = /(*) " /( i*)- 



r(a)( 2% — x ) 



1-1- 



The nonholonomic geometry of fractional tangent bundle depends on the 
type of chosen fractional derivative. We also emphasize that the above formu- 
las can be generalized for an arbitrary vector bundle E and/or nonholonomic 
manifold V. 
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A. 2 N— and d— connections and metrics 

a a 

A nonlinear connection (N-connection) N for a fractional space V is 
defined by a nonholonomic distribution (Whitney sum) with conventional h- 

a a 

and v-subspaces, hV and vV, 



a a 



TV = hV®vV. (A.4) 

a 

A fractional N-connection is defined by its local coefficients N={ a N^}, 
when 

N= a ATf( M )(^)«®§ a . 

a 

For a N-connection N, we can always construct a class of fractional (co) 
frames (N-adapted) linearly depending on a iV°, 



*e, = 5, - a Nfd a , a e b = d b 

a fS _ f a J _ (J„j\a a„b _ (J..b\a , a Ar b(j^k\a 



(A.5) 



e" = [ a e J = (dx J ) a , V = (dy b ) a + a N b k (dx k ) a \. (A.6) 

a 

The nontrivial nonholonomy coefficients are computed a W® b = <9 fe a Nf and 



°Wg = a fi^ = "e^ a A^ a - %• a A"f for 



In above formulas, the values a f^ are called the coefficients of N-connection 

curvature. A nonholonomic manifold defined by a structure N is called, in 
brief, a N-anholonomic fractional manifold. 

a a 

We introduce a metric structure g = { a gap} on V as a symmetric 
second rank tensor with coefficients determined locally with respect to a 
corresponding tensor product of fractional differentials, 

g= a g J i{u){dvX) a ®{duty. (A.7) 

For N-adapted constructions, it is important to use the property that any 

OL 

fractional metric g can be represented equivalently as a distinguished metric 
(d-metric), g = [ a g kj , a g cb ] , when 

g = a g kj (x,y) a e k ® V + a g ch (x,y) a e c ® a e b (A.8) 

= Vk'j' a e k ' <g> V + T) dV a e c ' <g> a e b ' , 
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where matrices r\ k iy = diag[±l, ±1, ..., ±1] and r) a iy = diag[±l,±l, ...,±1], 
for the signature of a "prime" spacetime V, are obtained by frame transforms 
Vk'f = e k k > eJ j> a 9kj and r) a >v = e a a , e b b , a g ab . 

a a 

A distinguished connection (d-connection) D on V is defined as a linear 
connection preserving under parallel transports the Whitney sum (1A.4J) . We 
can associate a N-adapted differential 1-form 

Q r r „= a T T ^ a e\ (A.9) 



parametrizing the coefficients (with respect to (1A.6J) and (1A.5J) ) in the form 

ap7 I aji aja ar^i a/^<a\ 

1 t8 — { ^jki ^bki °jc) °6cJ ■ 



The absolute fractional differential a d = lX d x + iy d y acts on fractional 

8 



differential forms in N-adapted form; the value Q d := "e 13 Q e« splits into 



exterior h- and v-derivatives when 

X X := (dxT ia J, = V "e,- and iy d y := (dy a ) a lX d a = V a e b . 

a 

The torsion and curvature of a fractional d-connection D = { a T T g } can 
be defined and computed, respectively, as fractional 2-forms, 

a T T = D a e T = a d a e T + a T T /3 A a e p and (A.10) 

a r 

i^ p — u i. p — a i ^ — i- p/\ l 7 — «- y3 7< j e a e. 

There are two another important geometric objects: the fractional Ricci 
tensor a 1Zic = { a ~R a/ 3 = Q R T Q , /3r } with components 

a p _j_ a pfc a p _i_ « pfc a p _j_ a pb a p _j_ a pc /A i i \ 

a 

and the scalar curvature of fractional d-connection D, 

^R= a g T/3Q R r/3 = "#+ a 5, "# = V J ' a Rij, a S = a g aba R ab , (A.12) 

with a g T/3 being the inverse coefficients to a d-metric (1A.8|) . 
We can introduce the Einstein tensor a Sns = { a G a/ 3}, 

a G a /3 '■= "R-o/3 — ~ "gad "R (A. 13) 

For various applications, we can considered more special classes of d- 
connections: 
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There is a unique canonical metric compatible fractional d-connection 
«D = { °T\ = ( "L* k , a L a bkJ °& j& Q Q C )}, when «D ( «g) = 0, 

satisfying the conditions a f ) k = and a f a bc = 0, but a f * ja , a f a . { and 

a T a bi are not zero. The N-adapted coefficients are explicitly determined 
by the coefficients of (1A.8J) . 

a T i — a n ir I a a a n _l_ a a a n a a a n \ 

L jk — 7) 9 { e k 9jr + Gj gkr ~ &r 9jk) i 

a^ra _ a. I a im\ , l a act a a 
L bk - e b{ ^k ) + 2 9 { e k 9bc 



9dc e b l\ k - g db e c l\ k ), 

1 

2 



a C) c = ^ a g ika e c a g jk , (A.14) 



uPia _ l a ad ( a a , a a a a \ 

<^bc — 2 9 { e c 9bd + e c g cd - e d gbc) ■ 

• The fractional Levi-Civita connection a V = { T 7 „} can be defined 
in standard from but for the fractional Caputo left derivatives acting 
on the coefficients of a fractional metric (IA.7I) . 

On spaces with nontrivial nonholonomic structure, it is preferred to work 
on V with a D = { T 7 T „} instead of a V (the last one is not adapted to the 
N-connection splitting ( 1A.4I) ). The torsion a T T ( 1A.10I) of a D is uniquely 



induced nonholonomically by off-diagonal coefficients of the d-metric ( 1A.8I) . 
With respect to N-adapted fractional bases (1A.5[) and (1A.6J) . the coeffi- 



cients of the fractional Levi-Civita and canonical d-connection satisfy the 
distorting relations 

^ = «f \ p + «Z\^ (A.15) 
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where the N-adapted coefficients of distortion tensor Z 1 B are computed 



ct/~ii a „ a. ab 



x ryi r\ a rva ot/~ii ot a ab a O a 

^jk — u > ^jk — ~ ^jb 9ik 9 - ^ "jfc, 
2 "Ik a 9cb a 9^-\w 



2 

am _ ' one a a ji x ( xi xh a a ft\ ar~ij 

Z bk - o ll jk 9cb 9 ~ 7,K d j d k ~ 9jk 9 ) <^ hb , 



a Kk = \(* a A+ a 9c d a 9 ab )[ a L c bk - a e b ( «N c k )}, 
a Zl b = \ a V% a g cb a g ji + l(6}6 h k - a g jk a g lh ) a C j hb , 



2 ° K J ^ J 2 
Z% = -\{8 a X- a g cb a g ad )[ a L%- a e d { a N<)] , *Z a hc = 0, 



a rya 



a Z ab = ~^{[ a L c aj - a e a {*Nf)] <*g cb 
+ [ a L c bj - "e b ( a N<)} a g ca }. 
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